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In a previous work [Phys. Rev. A 66, 0134XX (2002)] we noted a partial disagreement between 
quantum i?-matrix and semiclassical calculations of photoabsorption spectra of molecules in a mag- 
netic field. We show this disagreement is due to a non- vanishing contribution of processes which are 
forbidden according to the usual semiclassical formalism. Formulas to include these processes are 
obtained by using a refined stationary phase approximation. The resulting higher order in h con- 
tributions also account for previously unexplained "recurrences without closed-orbits". Quantum 
and semiclassical photoabsorption spectra for Rydberg atoms and molecules in a magnetic field are 
calculated and compared to assess the validity of the first-order forbidden orbit contributions. 

p • PACS numbers: 32.60.-l-i, 33.55.Be, 03.65.Sq, 05.45.-a 



. I. INTRODUCTION 

^ The photo-absorption spectrum of excited atoms or molecules placed in a magnetic field displays complex structures. 



Closed-orbit theory consists of a fully quantitative approach in which the large-scale structures of the spectra are 
explained in terms of classical trajectories closed at the nucleus, i.e., leaving and returning to the core. Each orbit 
produces on its return an oscillation in the photo-absorption cross-section; the Fourier-transform of the spectrum, 
known as the recurrence spectrum therefore exhibits sharp peaks at the period of the orbits. First developed for 
' the hydrogen atom [Q, closed-orbit theory was then extended to treat the case of non-hydrogenic Rydberg atoms: 
\ the additional ^ectral modulations appear as the result of successive quantum encounters of the Rydberg electron 
I with the core . The wave-function follows the hydrogenic classical orbits in the region where the Coulomb and the 
■ external fields compete ("outer region"), but near the core, where the external fields are negligible ("inner region"), 
[ — : the wave-function is described quantum-mechanically. 

' More recently, we proposed a closed-orbit theory treatment of molecules in external fields in addition to 

, elastic scattering, inelastic scattering gives rise to novel spectral modulations. The relative importance of elastic and 
inelastic scattering was shown to depend on the short-range phase-shifts, the molecular quantum defects. In the 
c/3 inelastic collision process, the molecular core undergoes a transition from its ground state to an excited state, and 
. ^ the dynamical regime of the Rydberg electron changes accordingly, say from a chaotic to a near integrable classical 
regime. Good quantitative agreements between quantum calculations and closed-orbit theory in the case of an external 
magnetic field were obtained. However we noted in Ref. [Q a discrepancy between the quantum and the semiclassical 
results for certain peaks of the recurrence spectrum. In general, disagreements between semiclassical and quantum 
recurrence spectra are due to higher order (in h) effects such as bifurcations or ghost orbits, and specific formulas to 
account for these effects have been developed |^ . The discrepancies we observed in are due to another type of 
effect, namely the manifestation of orbits which are forbidden to first order in h; these orbits are forbidden either i) 
5^ ■ because they should not be excited according to the usual semiclassical formalism or ii) because they do not classically 
' exist. 

We investigate in this work the effects of such first-order suppressed orbits. We will derive formulas for including 
their contribution in the recurrence spectra by going beyond the usual stationary phase approximation employed in 
the standard form of closed-orbit theory. The formulas will be tested versus exact quantum calculations for different 
model atomic and molecular systems in an external magnetic field. The two orbits that will be dealt with specifically 
are the orbit perpendicular to the field, which should not be excited when it lies in the node of a wavefunction, and 
the orbit parallel to the field, which does not exist classically when the electron's angular momentum projection on 
the field axis is non- vanishing (since the Hamiltonian remains divergent on the z axis even after regularization). 
The contribution of orbits lying in the node of a wave-function was first observed by Shaw et al. [Q when comparing 
quantum and semiclassical calculations for the diamagnetic hydrogen atom in the near-integrable regime (scaled 
energies e ~ —0.7). They obtained a formula for including the first-order forbidden contribution of the perpendicular 
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orbit, which appeared as a small feature in the recurrence spectra. In non-hydrogenic systems, we can expect the 
effects associated with "forbidden" orbits to be far more important than in hydrogen given that core-scattering mixes 
the contributions of different orbits. Moreover, although it may have been expected that at higher scaled energy 
the contribution of on-node orbits would become insignificant (as the classical amplitudes decrease), we will see that 
their inclusion is necessary to account for the correct amplitude in the modulations produced by orbits which have 
bifurcated from them. The existence of recurrences produced by "non-existing" orbits on the field axis was reported 
for Rydberg atoms in an electric field by Robicheaux and Shaw |8| , who developed a heuristic formula which yielded 
a poor agreement between semiclassical and quantum calculations. We will derive a formula for appropriately taking 
into account such classically forbidden orbits and compare it to quantum calculations in the case of an external 
magnetic field. In passing we will also show that the contribution of the parallel orbit when it is allowed (i.e., for 
= 0) can be obtained by treating it as any other orbit, provided a higher order refined stationary phase integration 
is used (whereas the parallel orbit has always been treated as a special case, following the original derivation given 
by Gao and Delos |^). 

The paper is organized as follows. We recall in Sec. II the usual semiclassical formulas of closed-orbit theory (with 
provision for multichannel core scattering). Sec. Ill details the derivation of the contribution to the photoabsorption 
spectra of the first-order forbidden orbits. We report in Sec. IV quantum and semiclassical calculations for different 
values of the quantum defects, scaled energies or magnetic field ranges, focusing on the contribution of those forbidden 
orbits. We give our conclusions in Sec. V. 

II. OVERVIEW OF STANDARD MULTICHANNEL CLOSED-ORBIT THEORY 

Closed-orbit theory explains the dynamics underlying the photoabsorption spectra of Rydberg atoms or molecules 
in external fields in terms of closed orbits: following initial photo-excitation, the wavefunction of the excited electron 
propagates first in a region near the ionic core ("inner region"), in which the external field can be neglected. Beyond 
the inner region, the wave- function is propagated semiclassically along classical trajectories. Some trajectories return 
to the inner region, and the semiclassical wavefunction carried by those trajectories is matched to an exact wave- 
function in the inner region given by a standard (field-free) multichannel quantum defect theory (MQDT) expansion. 
The superposition of these returning waves with the initially dipole-excited wavefunction produces sinusoidal modu- 
lations in the photoabsoprtion spectrum, which appear as isolated peaks in the Fourier-transformed ( "recurrence" ) 
spectrum. Further modulations (i.e., peaks in the recurrence spectrum) are caused by the core-scattering process; in 
a multichannel problem, the electron can exchange energy and angular momenta with the core, so after the collision 
the electron wavefunction propagates outward again, and when it leaves the inner region the wavefunction will follow 
once again classical trajectories. If the collision was perfectly elastic, the electron will follow one of the previously 
followed trajectories; on the other hand, inelastic collisions will result in trajectories pertaining to a different classical 
regime. 

In Ref. Q we described in detail photoabsorption from a ground state diatomic molecule in a static magnetic 
field. After photoexcitation, the molecular core could either be rotationally excited [N = 2) or non-excited {N = 0, 
where N is the core angular momentum). The molecular core then plays the role of an effective 2-level scatterer 
which combines classical trajectories belonging to 2 different dynamical regimes (typically, chaotic and near-integrable 
regimes), thereby producing additional modulations in the photoabsorption spectrum. The way these combinations 
occur depends both on the classical characteristics (amplitude and action Sk of the fcth trajectory), and on the 
properties of the scatterer (which is given by the scattering transition matrix, T) . Scaled energy spectroscopy consists 
of simultaneously varying the magnetic field strength 7 and the laser excitation frequency so as to keep e = E^~'^^^ 
constant, where E is the energy of the Rydberg electron; e is the scaled energy, which depends on the core state 
j through the energy partition between the core and the outer electron. Although scaling for molecules is only 
approximate, we have seen in Q how a molecular system can be scaled conveniently; will stand for 7^^'^, since the 
field strength plays the role of the Planck constant jlOj . The absorption rate in the one core-scatter approximation is 
then given by [see Eq. (3.30) in [||] 

Tih^s) = 21^/^3/^ ^^^Im{(a| j) C„C„,e'-(''°+'^=') 

j a a' 

{j\ a') Y: Kie,) + ;.^/^2"/4,3/2 J2 (/I a') T,, ^ n{ {e,) ^ UHe,) 

k j' k q 

The notation has been completely detailed in Sec. Ill of |4|, but in short: j (and j') is a compound index accounting 
for the core-electron couplings when the electron dynamics has uncoupled from the core. If we assume H2 as the 
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prototype molecule, the initial state has the quantum numbers J = 0, / = ( J is the total angular momentum, I the 
orbital momentum of the outer electron), the sum over j runs over the core states N = 0, Mj^ = M and N = 2, 
Mn = 7\/ — 1, A/, A/ + 1 where Mn is the projection of N on the field axis; M, the projection of the total angular 
momentum on the field axis, is the only quantum number conserved throughout the entire physical process, a gives 
the set of quantum numbers in the molecular frame, when the electron is coupled to the molecular axis; the sum runs 
here on A = (S state) and A = 1 (11 state); /is and /in are the corresponding short-range molecular quantum defects. 
Cq is a coefhcient giving the relative strength of the electronic dipole transition amplitudes; from the united dipole 
approximation, known to be valid for H2, we have Cs = 1 and Cn — a/S- {j\ 01) are the transformation coefficients 
between the molecular and the uncoupled frames. The elements of the scattering matrix Tjji depend solely upon 
the quantum defects and the (j| a) elements. The quantities 7?.^(ej) are the only ones that depend on the classical 
properties of the Rydberg electron trajectories. We have, for the trajectory k associated with the core in state j [Eqs. 
(3.18) and (Dl) of g]: 

^i(e,) = |sin0,fe sin0/fc|i/2 ^(-l)'^+'.'ll^.,™^., {e,k)Yi*^^ {Ofk) 



where 1^,1]' ^ |to, | and m,- 



9ik and 9fk are the initial and final angles of the fcth trajectory relative to the 
magnetic field direction, which is taken to be along the z axis. A;'^"^ and ^ ^ are the scaled classical amplitude 

and action, evaluated at the corresponding scaled energy e^; uj^^™^ is the associated Maslov index. Yi-mj {Oik) will be 
used throughout as a short-hand notation for 11 m (^ifc,0), since the conserved axial symmetry has been separated 
from the 2-dimensional scmiclassical problem; we have accordingly used a different notation for the quantized value 
of the Rydberg electron's angular momentum projection m and its classical counterpart Lz appearing in the two- 
dimensional diamagnetic Hamiltonian. Equation ( 2|) must be modified for orbits lying along the magnetic field axis 
{Oik = Ofk = 0) as detailed below. Note that Eq. (n ) is also valid for ground state hydrogen photoexcited to odd-parity 
states (by setting /is and /in to zero; then the T matrix vanishes) as well as for non-hydrogenic Rydberg atoms with 
a single quantum defect /i;=i (by setting /is = /in — Mi=i; the T matrix is then diagonal). 

Eqs. (|l|) and (||) are obtained by matching the scmiclassical wavefunction V'gc"™^ associated with the core in state 
\Njmj) which returns to the core region to a MQDT expansion on a boundary circle {rf,6f). The scmiclassical 
wavefunction reads 



i^^s'P{rf.9s) 



J2^or{r: 

k 



i, Oik) 



rf sin 0ik 



rj sm U f 



1/2 



{rj,ej) expz (^f^"^ {rj^ej) - cf^"^ n/2) 



(3) 



where ipout^^ (''ii ^ik) represents the initially outgoing waves, which were propagated semiclassically beyond the bound- 

, AT, 



ary (r,, 6*,). We shall write V'out ' {^i, ^ik) as 



(4) 



lj^\mj\ 



with 



(5) 



where the are the dipole transition amplitudes in the molecular frame. The MQDT expansion reads in the 
uncoupled basis 



(6) 



where / and (/"*" are Coulomb functions (/ is regular at the origin, is an outgoing wave). The expansion coefficients 
Cji are obtained by matching Eqs. (||) and (^ on the boundary. The matching condition reads 



(r/)/2i = -27r / d0f sinOfY^* {9f)t^^^"^^ {rf,0f) 



(7) 
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the integral is performed for each traject ory k in the stationary phase approximation^ since the phase is stationary 
along the final angle of the trajectory 9fk The value of the coefficients Cj are then inserted in the expression giving 
the dipole transition amplitudes, of the form (V'ol D |'(/'qdt) where \4>o) is the initial state prior to photoabsorption. 
The formulas for the oscillator strength and the absorption rate are then obtained. Obviously when Eqs. (||) or (|^) 
vanish, e.g., Yi.raj{Oik) = or Y-i*^.{0f) — 0, then TZ'^ej) vanishes, and these orbits should not produce modulations 
in the oscillator strength. 



III. CONTRIBUTION OF FIRST-ORDER FORBIDDEN ORBITS 



A. General remarks 

Although Eqs. (P and predict that if 'Rjj^[e.j) vanishes, the orbit k should not contribute to the recurrence 
spectrum, we had observed in 1 1| a mismatch between semiclassical and quantum classical calculations for molecules 
in fields in the amplitude of certain peaks in the recurrence spectrum. This mismatch was interpreted as arising 
from the interference of the orbit perpendicular to the field (which lies on the node of a spherical harmonic when 
M = and is thus semiclassically forbidden) with the R\ "pac-man" orbit. As stated, the recurrences associated with 
classical orbits lying in the node of a wavefunction were first observed by Shaw et al. when comparing quantum 
and semiclassical calculations for the diamagnetic hydrogen atom at low scaled energies (e ~ —0.7). They obtained a 
formula for the contribution of the perpendicular orbit by matching the returning semiclassical wave to an "ansatz" (a 
rotated first order Bessel function) . In this section we shall derive simply the contribution to the oscillator strength of 
this type of first-order suppressed orbit by employing the same framework introduced in Sec. II, without introducing 
additional assumptions; only the stationary phase integration needs to be performed differently. We will also derive a 
formula to account for the peaks in the m = 1 recurrence spectra appearing at the scaled action of the m — Q parallel 
orbit, since the parallel orbit does not exist classically when Lz ^ and there is therefore no corresponding TZj,{ej) 
factor. Numerical results and examples will be given in Sec. IV. 



B. Contribution of on- node suppressed orbits 



The rationale for including the contribution of orbits lying on the node of a wavefunction was already given in [Q: 
strictly speaking, an orbit k closed at the core with initial and returning angles 9ik and 6fk is not isolated, but has 
neighboring orbits which are not closed at the origin. We assume a neighboring orbit returns with an angle 9f and 
envisage the initial angle di of this orbit to be a function of 9f, i.e., 6i = 0{9f). To first order in 9f — 9fk we have 



Yl.mj{Qi) — Yi.mj{9ik) + 



dYi^^^{9,k) d9,k 
d9,k 



89 i 



'Ik 



89 



f 



(8) 



(9) 



In the usual case, the contribution of the neighboring orbits with initial and final angles {9i,9f) is negligible when 
compared to the central orbit with angles {9ik,9fk)- However, when the central orbit lies on a node of a spherical 
harmonic, the semiclassical wave-function can only be carried by the neighboring orbits, and this contribution can be 
significant provided the classical density of trajectories is sufficiently large on return to the core. 

It turns out that, to first order in 9f ~ 9 fk, such a contribution comes into play if we have both Yi^rnj {9ik) ~ and 

Yim- i^fk) — 0- If ^ is such an orbit, its contribution to the outgoing wave ipout^^ [E'l- ®l] vanishes. The contribution 
of the neighboring orbits are taken into account by inserting Eq. (^) in Eq. (||) and Eq. (||) in Eq. (^; the right 
hand-side of Eq. (j^) then takes the form: 



-2tt 



d9f |sin6'f sin 9ik 



1/2 



E2V 



^Yl^,rr.,{^^k) 39, k 



89,k 



39 f 



dYi*^,iOf) 



39 f 



'fk 



(10) 



Following Hiipper et al. ||T^, we express the action on the boundary (rf,9f) in terms of the action of the orbit closed 
at the origin, S^^™^ {rf,9f) ~ + \/''//8 {9f — 9fk) ■ The integral can now be performed; a straightforward 
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stationary phase integration would lead to zero, since the integrand vanishes at the point of stationary phase 9f = 9fk. 
However, we can assume the integrand to vary slowly around the angle of stationary phase, and integrate exactly the 
term between square brackets (see Appendix A). In the semiclassical limit, Eq. (A4) is appropriate. Eq. ( |To| ) then 
becomes 

d9a: . 



fc i7r/2r,l/2 — l/4<A''jm3/ ^ s / 'J"ik \ 

fte"'^2^'^ A,,' {rf,9fk) sgn(^^) 

Isin^., sine,,^/^ rJ^/^Af™^ (r,, 0,,) exp [^ {s^^i^,^ ~ 2^ - -f'"^V2)] } 



d9fk 



(11) 



where we have used A^^™^ {rf,9fk) = ■ For clarity we have singled out the factor specific to on- node orbits 

(in front of the curly brackets) relative to the expression valid for "typical" allowed orbits (inside the curly brackets) . 
In particular it can be seen on-node orbits are suppressed by a factor H relative to typical orbits. 

The relevant scaled factor ^fe_nodc(^j) S^^i^^S contribution of an orbit lying on the node of a wavefunction in 
the absorption rate is thus 



1/2 



■'^fc-nodcV 



(e,) = ?icff2i/2sgn(^)|sin0,fcsin%|i/'^(-l) 



'd9fk 



d9ik 



d9f 



exp I 



(12) 



This formula holds for non-vanishing angles. Specializing to our atomic and molecular model described above, we 
have / = 1, so this formula only applies to the orbit perpendicular to the field (^jj. ~ 9fk — 7r/2) for the m = 
manifolds (e.g., for a molecule, when M — 0, for an outer electron associated with core states having a projection 
Mjv = 0). Note that in the absence of core-effects, Eq. (|l^) becomes strictly equivalent to the correction obtained in 
for the hydrogen atom. 



C. Contribution of the "classically non-existing" parallel orbit 

1. Contribution of the parallel orbit when Lz = 

We first recall that the orbit parallel to the field {9ik — 9fk = 0) classically exists if = 0. Even then, this orbit is 
treated as a special case, because the formulas valid for the other orbits, Eqs. (0) and (||), need to be modified. This 
modification was originally obtained by matching the semiclassical returning wave to a particular Bessel function on 
the z axis ^ . We show here that the reason this orbit is "special" is that the standard stationary phase approximation 
vanishes. Indeed, setting \shi9ik/ sm9 fk\^^^ — > \d9iii/d9 fk\^^^ in the outgoing wave (|^), the expression to be integrated 
arising from the matching condition Eq. (Q) is 

^ d9 f sin 9 fe:>q){i^lVJJd,9}/h), (13) 

which is zero in the standard stationary phase approximation. However, an approximate closed form may be obtained 
(see Appendix B). To first order in h, we have 

^ d9fsm9fexp{i^jrf/89j/h) « He'^'/^y^. (14) 

This result is reduced by a factor /i^/^e'"^/^2~^/''7r~^/^ relative to the standard stationary phase integration for non- 
zero degree orbits, which is exactly the result obtained in . The parallel orbit thus appears as a first-order suppressed 
orbit, which is apparent from its /i-dependence. 
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FIG. 1: Recurrence spectrum (Fourier transform of the photoabsorption spectrum) for a non-hydrogenic atom with = 0.5, 
M = 0, at e = —0.7, in the range — [60, 120]. Top: quantum calculations. Bottom: standard semiclassical calculations 

(solid line), semiclassical calculation including the higher-order contribution from the on-node orbit (broken line). 



2. Contribution of the parallel orbit when ^ 



When Lz is non- vanishing the diamagnetic Hamiltonian contains the repulsive term proportional to Ll/fP, where 
Lz = j^^^Lz is the scaled angular momentum and p = j^^^p the scaled distance from the z axis 11^]. Even though 



in the semiclassical limit Lz is small (since hcs = 7 



1/3 



0), the centrifugal term is infinite on the z axis, and the 



parallel orbit no longer exists |1^. However, we may expect orbits neighboring the z axis and not closed at the 
nucleus to contribute to the oscillator strength, in the same manner as for the on-node orbit (orbits near the z axis 
for Lz ^ and their structural stability as — > where actually investigated in |l2j). Starting from Eqs. (||) and 
(^ and setting |sin6'ife/ sin6'/fe|"'^^^ \d9ik/d9fk\^^^ in Eq. (||) as in Sec. Ill B leads, after matching the semiclassical 
returning wave to the MQDT expansion, to the integral 



^6*/ sin 6*/ 6*; exp (i J r/ /89j /K) 



To lowest order in /i, we have (see Appendix B) 



dO / sin 6*/ 0j exp{iJrf /Wj /K) 



/o ' r 

The resulting scaled contribution to the oscillator strength is given by 



(15) 



(16) 



dOik 
exp i 



(17) 



Within our molecular model, this correction applies when |m| = 1; this is of course the case when Af = 1, but 
even for M = 0, the outer electron may be associated with core states having a projection |M7v| — 1, i.e., N ^ 2 
Mjv = 1,— 1- Note that we have followed the by now standard notation whereby the zero-degree orbit amplitude is 
'fk = 0) = \d9ik/d9fk\ (although stricto sensu this is the square of the genuine two-dimensional 

rf,dfk = 0) is independent of the boundary radius r/. 



set as A^'™'^{rf 

semiclassical amplitude), so that now f^.^^^^f''™' 



D. h dependence 



Unsurprisingly, the contribution of the forbidden orbits in the recurrence spectra have a different h dependence. 
The on-node orbit is suppressed by a factor h relative to a typical primitive orbit; the parallel orbit is suppressed by a 
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FIG. 2: General view of the recurrence spectrum for a molecule with the set of quantum defects fis ~ —0.3, fin = 0.1, M — 0, 
at £jv=o = —0.3 and ejv=2 = —0.8 in the range 7"^''^ = [60, 120]. The semiclassical result (bottom) has been calculated in the 
one core-scatter approximation, but the first-order suppressed contributions have been included. 



factor h}^'^ relative to a typical orbit, and the forbidden parallel by a factor h relative to the classically allowed parallel 
orbit and h'^^'^ relative to a typical orbit. This is to be contrasted with the core-scattered ( "diffractive" ) orbits: each 
encounter with the core brings in for a typical orbit a factor h^^"^ . Thus single core-scattering is expected to dominate 
the photoabsorption spectrum in the semiclassical regime; however, the h dependence is balanced by the amplitude 
factors, explaining why for individual orbits the forbidden contribution may be strong, as will be seen below. It 
may also be noted that the combination of orbits having different individual h dependence through core-scattering 
[last term in Eq. (^] will give rise to peaks in the recurrence spectra with a dependence of the form h'^^^, where 
V is an integer depending on the type of primitive orbits connected by the core-scattering process. In particular, 
core-scattering between two forbidden parallel orbits is expected to be highly suppressed in the semiclassical limit. 



IV. RESULTS 



We compare below quantum and semiclassical calculations to assess the importance of the forbidden orbits in the 
recurrence spectra of atoms and molecules. The numerical examples given in this section correspond to non-hydrogenic 
atoms and different molecules obtained by choosing different sets of quantum defects, within the framework of the 
model described in Sec. II. 

Fig. 1 displays the recurrence spectrum of a non-hydrogenic atom with — 0.5, M — 0, at e = —0.7, in the range 
7"-'^/^ = [60, 120]. The top figure gives the quantum calculation, whereas the solid line in the bottom part of the plot 
results from the standard semiclassical treatment; this solid line only accounts for less than half of the peaks in the 
recurrence spectrum. The missing peaks relative to the quantum results arise from the orbit Ri perpendicular to the 
field (and its nth repetition i?„) — which lies on the node of the wavefunction and is thus not excited according to 
the standard treatment — as well as from the combinations produced by core-scattering between i?„ and the parallel 
orbit and between the on-node orbits. The broken line includes the contribution of the on-node orbit [Eq. (|l^)] in 
the semiclassical calculation. 

At higher scaled energies, the contribution of the forbidden i?„ orbit is visible through the mismatch observed in 
[Ql between the height of the peaks in the quantum and semiclassical recurrence spectra. Fig. 2 displays a global 
view of the recurrence spectrum for a molecule with the set of quantum defects /is = —0.3, /in = 0.1, M = 0, at 
iN=o — —0.3 and €n=2 = —0.8 in the range 7~i/-^ = [60, 120]. These quantum defects yield a balanced contribution of 
the different type of orbits: the primitive geometric orbits (that is the orbits that appear in the recurrence spectrum 
of the hydrogen atom), the elastic scattered diffractive orbits (that appear in the recurrence spectra of non-hydrogenic 
atoms and in molecules) and the inelastic scattered diffractive orbits (that solely appear in molecular systems). In 
Figs. 3-5, we zoom on some individual peaks in the recurrence spectra, choosing different sets of quantum defects but 
keeping the other parameters (scaled energies, 7 range) constant, to observe the presence of the on-node orbit and 
how its interplay with core-scattering affects the amplitude of the recurrence peaks. 

Fig. 3 displays the recurrence spectra for non-hydrogenic atoms with = 0.5 [a], fj-i^i — 0.25 [b] and fj.i=i = 0.1 
[c], at e = —0.3, around the peaks labelled 3 and 1 + 3 in Fig. 2. According to the standard treatment (solid line), 
peak 3 is produced by the i?2 "pac-man" orbit (the shapes and characteristics of the orbits mentioned here are given 
in Table I and Fig. 6 of [Q ; i?2 has bifurcated from R2 at a slightly lower energy, and thus the two orbits have nearly 























CO 










un 


















b / 






E 






X 


o 








iransl 








1 








'ouri 






X 











3 



3 



1.95 3.9 

Scaled Action 



FIG. 3: Recurrence spectra for non-hydrogenic atoms with fii^i = 0.5 [a], = 0.25 [b] and /i;=i =0.1 [c]. The left panel 
represents the peak labelled 3 in Fig. 2, the right panel shows the combination peak 1 + 3 (the amplitude of the 1 + 3 peak has 
been multiplied by 3 relative to the amplitude of the peak on the left panel) . For each peak, the quantum result (top) is plotted 
versus semiclassical calculations (upside-down) without (solid line) and with (broken line) the higher order contributions. 

the same scaled action), and 1 + 3 results from the combination of the "balloon" orbit (peak 1) and i?2 through 
core-scattering. The mismatch for the peak 3 arises from interference between the contributions of the i?2 and the 
on-node i?2 orbit; indeed, including the on- node orbit in the semiclassical calculations results in excellent agreement 
with the quantum result. The peak 3 thus results from the interference of primitive orbits and accordingly does not 
depend on the value of the quantum defect; however the peak 1 + 3 does depend on the quantum defect and vanishes 
in the limit — *■ 0; the contribution of the on-node orbit in 1 + 3 is seen to be important (in absolute terms) only 
provided the quantum defect is large. Note that in principle we should also have taken into account the first and 
third repetitions of the perpendicular orbit, but their corresponding amplitudes are very small, so these orbits have 
a negligible contribution to the recurrence spectra. 

The situation depicted in Fig. 4 is more involved: a close up of the peak at S' = 3.9 (labelled 3 + 3 in Fig. 2) 
is shown for a molecule with quantum defects /zs = 0.5, /xn = [a] and = 0.5, = 0.5 [b]; the peak 3 + 3 
arises from recurrences produced by different orbits: the second return of i?2 and the fourth return R4 of the on-node 
perpendicular orbit, the combinations R^ + R2, R2 + R2 and R2 + R2 via core-scattering. The resulting peak amplitude 
depends both on the quantum defects (which rule the core-scattering amplitudes) and on the inclusion of the two on- 
node orbits: in the first case the standard semiclassical result underestimates the exact quantum calculation, whereas 
in Fig. 4 (b) the standard semiclassical result overestimates the correct recurrence strength. Adding the contribution 
of the on-node orbits in the semiclassical treatment results in both cases in a better agreement with the quantum 
calculations. 

Fig. 5 displays the peak labelled a + 3 in Fig. 2 but for the choice of quantum defects /zs = 0.5, /in — 0. This peak 
results from the inelastic scattering between R^ aX e ^ —0.3 and the perpendicular orbit associated with the core 
state N = 2, m = ±1 at e = —0.8. Again, the standard closed-orbit result underestimates the recurrence strength 
and the inclusion of the first-order suppressed on-node orbit improves the agreement with the quantum results. 

Finally, Fig. 6 shows a portion of the recurrence spectrum for the hydrogen atom at e = —0.55, M = 1, in the 
range 7"^/'^ = [30, 240]. We have zoomed the peaks at = 0.95 and S = 1.91 which are due to the first and second 
returns of the classically "non-existing" forbidden parallel orbit. Note that the peak at 5' = 0.95 sits on the right 
shoulder of the much stronger Ri orbit, whereas the second return at S* = 1.91 is sufficiently isolated. The quantum 
calculation for M = 1 thus displays peaks for orbits which classically "do not exist" , at the actions of the corresponding 
M = parallel orbit. The standard semiclassical treatment (solid line) can not obviously account for those peaks, 
but including Eq. (p^), which takes into account higher order contributions, yields an excellent agreement with the 
quantum results, since those forbidden orbits contribute, albeit modestly, to the photoabsorption spectrum. 
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FIG. 4: The peak labelled 3 + 3 in Fig. 2 is shown for a molecule with quantum defects /is ~ 0.5, ~ [a] and /is = 0.5, 
fin — 0.5 [b]. The inclusion of the higher order contributions (broken line) gives a better agreement with the quantum 
calculations (top) than the standard semiclassical f^^v^-^i;^"^ ("o^i;^ 
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FIG. 5: The peak labelled a + 3, due to inelastic core-scattering in Fig. 2 is shown for the choice of quantum defects /ie = 0.5, 
fin = 0. The inclusion of the higher order contributions (broken line) gives a better agreement with the quantum calculations 
(top) than the standard semiclassical formalism (solid line upside-down). 



V. DISCUSSION AND CONCLUSION 



The feature developed in this paper is one of the many refinements that can be undertaken to improve a semiclassical 
formalism such as closed orbit theory. Some processes are forbidden on purely classical grounds (e.g., the above-barrier 
reflection of excited lithium in an electric field which results in very broad resonances in the absorption spectrum [ p^ ) 
whereas other processes are semiclassically distorted (e.g., diverging amplitudes at bifurcations). 

The role of the perpendicular on-node orbit was first observed in calculations for the diamagnetic hydrogen atom 
at low scaled energy (e ~ —0.7) 0. Subsequent high-resolution experiments on helium in a magnetic field in the 
same dynamical regime did not clearly detect the on-node orbits (they were within the experimental noise) 0. We 
have given a simpler derivation of the contribution of these first-order suppressed orbits, and our numerical results 
indicate that on-node orbits are more likely to be detected in non-hydrogenic atomic or molecular systems with strong 
quantum defects. At low scaled energies, peaks resulting from the core-scattering of the on-node orbit with a strong 
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FIG. 6: Recurrence spectrum for the hydrogen atom at e — —0.55, AI = 1, in the range 7"^^"^ = [30, 240]. The plot focuses on 
the first and second repetitions of the first-order suppressed "parallel" orbit, which is clearly visible on the quantum calculations 
(top). This feature is absent from the standard semiclassical calculations (bottom, solid line), but the inclusion of the higher 
order contributions (bottom, broken line) results in an excellent agreement with the quantum calculations. Note that the 
peak at S — 0.95 sits on the right shoulder of the much stronger Ri orbit; the oscillations are due to the finite range of the 
Welch-windowed Fourier transform. 



allowed orbit could be more easily detected; at higher scaled energies, the on-node orbit is most likely to affect the 
amplitude of peaks due to typical allowed orbits. 

The presence of contributions in the quantum photoabsorption spectra which were not correlated with any classical 
orbit was observed in calculations for non-hydrogenic atoms with m 7^ in an electric field by Robicheaux and Shaw 
; these contributions were coined "recurrences without closed orbits" because they appear at the scaled action of 
the parallel orbit which only exists classically when — and should therefore be absent in an m ^ recurrence 
spectrum. These authors also gave an ad-hoc semiclassical formula akin to the on-node correction which resulted in 
a poor agreement with the quantum calculations. Main p5| later pointed out that, for small but nonvanishing L^, 
periodic orbits having nearly the same action as the Lz = parallel orbit do exist; it was unclear however whether 
the "recurrences without closed orbits" could be attributed to such orbits, in particular because the starting point 
of these orbits is several atomic units away from the core. Our formula Eq. (|l^) correctly accounts for the peaks in 
the recurrence spectrum associated with these apparently non-existing orbits; the H dependence is different to that 
of the suppressed on-node orbits. The physical picture is similar in both cases: just as Eq. ( p^ accounts for close 
neighbors to the on-node orbit, which arc not closed at the origin but carry a portion of the wavcfunction back to 
the core region, Eq. (|l^) takes into account non-radial orbits close to the z axis which also give rise to recurrences by 
carrying the wavcfunction from and into the core region. 

To conclude, we have seen that first-order forbidden processes can be included within Closed-orbit theory in a 
simple and unified manner by elementary manipulations of the stationary phase integral, which yield a higher-order 
h dependence. In passing, we have shown that the zero degree orbit, which has always required special treatment, 
is in fact a case calling for a refined stationary phase integration. Analogous manipulations of the stationary phase 
integral of the Green's function were performed in to obtain an improved semiclassical long-range scattering 
matrix for Rydberg atoms in fields. Our method provides a convenient and effective way of including non-radial and 
non-closed trajectories that nevertheless contribute to the photoabsorption spectra of Rydberg atoms and molecules 
in fields without the need to to calculate explicitly the involved classical dynamics of those trajectories. The validity 
of the method was assessed by comparing our semiclassical results to quantum calculations for Rydberg atoms and 
molecules in an external magnetic field. 
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APPENDIX A 

We briefly work out the integral needed to determine the contribution of an orbit lying on the node of a wavefunction 
in Sec. Ill, 



(Al) 



This integral can be integrated directly but, for present purposes, it is convenient to express it in terms of sine and 
cosine Fresnel integrals and take the limit for the range in which the standard stationary phase approximation holds 
for the usual orbits. For example, in the neighborhood of 9f}~, the real part of (Al) can be expressed in the form 



I{n) ^ 2 



>/k+e„ 



_23/4ft3/2^-3/4^1/2 



{25(\/l + 2n) - 2Vl + 2ncos7rn} , 



(A2) 
(A3) 



where S is the sine Fresnel integral and £„ = 2^/'*r~^/^7r^/^?i^/^ vT~F2n with n a real number n > —1/2. For large 
half integer values of n, I{n + 1) — I{n) « and S{y/\ + 2n) ^ 1/2. I{n) can then be approximated by the term 
between square brackets in Eq. (AS). This is consistent with having neglected terms of order {0 f — 6 fkY Eq. ( Al) 
provided ?i — > 0. The imaginary part of Eq. (Al) is treated in the same way by writing the result in terms of the 
cosine Fresnel integral C{x). Hence 



defi0f - 0fkY eM^\ rflKQf - OfkY/K) 



Note that this result is independent of the value of 6'/^, provided 9fk 7^ 0. 



APPENDIX B 



j.3/42-5/4 ■ 



(A4) 



A closed form expression for the integrals 

h{e) 

and 



dO f sin Of exp(is0i 



He)= defsinef0jexp{isej), 



(Bl) 



(B2) 



with s real, are obtained in the limit in which the standard statio nar y phase approximation holds for the usual orbits 
by replacing the upper bound by e — > oo. Then Eqs. (Bl) and ( |B2| ) are given in terms of infinite series JTtI , which 
are actually representations of special functions. Choosing for simplicity a representation in terms of Fresnel integrals, 
Eq. (Bl) becomes in this limit 



/i(c») 



+ iS 



whereas for Eq. (B2) we have 



1 



^2(oo) = --3 1+ - (2s-z)exp 



1/2 



c 



(27rs) 



-1/2 



(27rs) 



+ iS 



-1/2 



(27rs) 



-1/2 



(B3) 



(B4) 



When s oo, to first order only C [(27rs) ^/^] ^ {2-ks) contributes to /i(oo) whereas for 12(00) the term between 
the braces simply gives 2. 
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